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ABSTRACT 

C"| , We present a catalogue of galaxy photometric redshifts and k-corrections for the Sloan 

Q_C Digital Sky Survey Seven Data Release (SDSS-DR7), available on the World Wide Web. 

The photometric redshifts were estimated with an artificial neural network using five ugriz 
bands, concentration indices and Petrosian radii in the g and r bands. We have explored 
^ 1 our redshift estimates with different training set concluding that the best choice to improve 

| redshift accuracy comprises the Main Galaxies Sample (MGS), the Luminous Red Galaxies, 

and galaxies of active galactic nuclei covering the redshift range < z ^ 0.3. For the MGS, 
the photometric redshift estimates agree with the spectroscopic values within rms = 0.0227. 
The derived distribution of photometric redshifts in the range < z p hot ^ 0.6 agrees well 
with the model predictions. 

' k-corrections were derived by calibration of the k-correct_v4 .2 code results for the 

, MGS with the reference frame (z = 0.1) (g — r) colours. We adopt a linear dependence of k 

(f) ■ corrections on redshift and (g — r) colours that provide suitable distributions of luminosity 

and colours for galaxies up to redshift z pno t = 0.6 comparable to the results in the literature. 
Thus, our k-correction estimate procedure is a powerful, low computational time algorithm 
capable of reproducing suitable results that can be used for testing galaxy properties at 
intermediate redshifts using the large SDSS database. 
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1 INTRODUCTION 

The knowledge of distances to galaxies is important to deduce 
intrinsic galaxy properties (absolute magnitude, size, etc.) from 
the observed properties (colours, sizes, angles, flux, apparent size, 
etc.). 

Since pioneering works measuring reds hift of bright galaxies (e.g., 
dShaplevI 4l933> . iHumason et alT l ll956h ). many efforts have been 
invested in mapping the light and matter in the Universe. Statistical 
analysis of galaxy properties and systems can be invaluable tools to 
study large-scale structure and evolution in the universe. 
In recent years, multi-band photometry has been performed for 
several millions of galaxies, whereas spectroscopic redshifts have 
been measured only for a small fraction of the photometric data. 
The Sloan Digital Sky survey has obtained multi-band images for 
approximately one hundred billion galaxies (and the next surveys 
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foresee increasing the number of objects up to billions), while spec- 
troscopic measurements have been obtained for nearly one million 
galaxies. A solution to the difficulty of obtaining spectroscopic red- 
shifts relies on the use of photometric redshift techniques. Although 
the redshifts calculated through these techniques are far less accu- 
rate than the spectroscopic measurements, these approximate dis- 
tance estimates allow for useful analysis in fields such as extra- 
galactic astronomy and observational cosmology. 

Two basic family of methods are commonly employed to 
calculate photometric redshifts. In the template matching approach, 
a set of spectral energy distribution (SED) templates is fitted to 
the observations (e.g., colours). In the empirical approach, on the 
other side, photometric redshifts are obtained from a large and 
representative training set of galaxies with both photometry and 
precise redshift estimations. The advantage of the first method is 
that it can also provide additional information, like the spectral 
type, k-corrections and absolute magnitudes. The accuracy of these 
estimations is limited by the SED models. The empirical model 
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overcomes this limitation through the use of a training set which, 
however, should be large and representative enough to provide 
accurate redshift estimations. 

The observed spectral energy distribution of distant galaxies 
is redshifted with respect to that in the g alaxy rest frame. The 
k-correction term dOke & Sandage 1 1 19681 ; iHogg et al.l 1200 ll) ap- 
plied to the apparent magnitude measured in a given photometric 
band takes into account this effect, allowing to compare photo- 
metric properties of galaxies at different redshifts. The estimation 
of k-corrections, then, is a requirement for many studies of dis- 
tant galaxies. It is possible to model k-corr ections as a function 
of redshift and galaxy morpholog i cal ty pe (Fukugita et alj 11995) : 
iMannucci et alj200lh . lLahav et alj ( fl995h and lBanerii et alJl201Cf) 
use ANNs to obtain morphological classification of galaxies. How- 
ever, these techniques employ as training sets, objects that had been 
classified by human eye subject to some degree of ambiguity and 
uncertainty, particularly at large redshifts. 

A more direct way to obtain k-corrections is by modelling 
galaxy SEDs as a function of wavelength. Usually template fitting 
of observe d galaxy fluxes is em p loyed to reconstruc t the S EP of 
the galaxy dBlanton et all d2003l) . lBlanton & Roweisl d2007h ). The 
feasibility and accuracy of this method relies in the quality of the 
models. 

For objec ts with spectral dat a, k-corrections can also be ob- 
tained directly. iRoche et al.l d2009l) used this technique to calculate 
k-corrections for early-type galaxies from the SDSS-DR6, provid- 
ing individual estimates for each galaxy. However, this technique is 

restricted to a limited n umber of galaxies w it h spec tro scopy. 

Recent works dChilingarian et alH d2010h , IWestra et al.l 

(2010)) have approximated k-corrections with analytical functions 
of re dshift, parametrized with s ome property characterizing galaxy 
type. Ichilingaria n et al.l d20ld) used different observed colour in- 
dices to approximat e k-corrections for nine filters {ugrizY J H K). 
IWestra et alj d2010h used spectra from the Smithsonian Hectospec 
Lensing Survey to obtain direct measurements of k-corrections by 
parametrization with the ratio of the average flux red and bluewards 
of the 4000A break (D n 4000). These kinds of parametrization 
simplify the computation of k-corrections and improve their 
accuracy. 

In this paper, we present a galaxy photometric redshift (z p hot) 
catalogue and a method for calculating k-corrections, for the sev- 
enth Data Release (PR7) of the Sloan D i gital Sky Survey ( SPSS ) 
imaging ca t alogue dBlanton et al. I 20031) , lEisenstein et al.l d200lh . 



York et alj d200d) ). To 



iGunn et alj dl998l) . IStrauss et al] d2002 
com pute photometric redshifts we used the ANNz software pack- 
age dCollister & Lahavl d2004h ). which have been shown to be 
a reliable tool. There are also two sets of photometric redshifts 
in the SPSS database: Abazaii an et al. I d2009h employ empirical, 
template-based and hybrid-te chniques app r oache s to photometric 
redshift estimation, whereas lOvaizu et alj [2008) adopt a neural 
network method and provide two different estimations, PI and 
CC2. Here we also com pare our redsh i ft est imations with those in 
the CC2 catalogue from lOvaizu et alj d2008h . which uses colours 
and concentration indices to infer redshifts. For the computation of 
k-corrections we propose a joint parametrization in terms of red- 
shift and the (g — r) colour in a certain reference frame for all 
SPSS bands, as well as an algorithm to determine them from the 
photometric data. We compare our results with those found in liter- 
ature. 

This paper is organized as follows. In Section 2 we describe 



the data used in our analysis. Section 3 presents our approach to cal- 
culate photometric redshifts in SPSS-PR7, analysing their advan- 
tages and limitations. Our estimation of k-corrections is presented 
in Section 4. Finally, Section 5 summarizes the results obtained in 
this work. 

Throughout this paper, we adopt a cosmological model char- 
acterized by the parameters Q, m = 0.3, SIa = 0.7 and Ho = 
75 ftkms" 1 Mpc" 1 . 



2 THE GALAXY SAMPLES 

The samples of galaxies used in this work were drawn from 
the Sloan Pigital Sky Survey S even Pata Relea se (SPSS-PR7, 
lAbazaiianetaLl d2009h ). SPSS dYork et alj d200d) ) mapped more 
than one-quarter of the entire sky, performing photometry and spec- 
troscopy for galaxies, quasars and stars. SPSS-PR7 is the seventh 
major data release, corresponding to the completion of the survey 
SPSS-II. It comprises 11.663 sq. deg. of imaging data, with an 
increment of ~ 2000 sq. deg., over the previous data release, ly- 
ing in regions of low Galactic latitude. SPSS-PR7 provides imag- 
ing data for 357 million distinct objects in five bands, ugriz, as 
well as spectroscopy over ~ n steradians in the North Galactic 
cap and 250 square degrees in the South Galactic cap. The aver- 
age wavelengths corresponding t o the five broad bands are 3551, 
4686 6165, 7481 and 8931 A dFukugita et all 19961 : IHogg et alj 
1200 it ISmith et alj|2002l) . For details regarding the SP SS camera 



Gunn et al. 1 1998); for astrometric calibrations see IPier et alj 

d2003l) . The survey has completed spectroscopy over 9380 sq. deg.; 
the spectroscopy is now complete over a large contiguous area of 
the Northern Galactic Cap, closing the gap that was present in pre- 
vious data releases. 

In this work we have extracted two galaxy data sets from 
SPSS-PR7, one with spectroscopic redshifts measurements and 
the other consisting in photometric data. 

The spectroscopic data was derived from the fits files 
at the SPSS home page; . consisting in Main Galaxy Sample 
(MGS: IStrauss et all 1 2002)), the Luminous Red Galaxy sample 



(LRG; Eisenstein et al. (2001)) and active galactic nuclei (AGN; 
iKauffmann et alj d2003h ). These spectroscopic samples were used 
for the computation of photometric redshifts and calibration of k- 
corrections. 

We built a first sample (hereafter Szl) consisting in 80% of 
the objects from MGS, 10% from LRG and 10% from AGNs, com- 
bined into a single set comprising ~ 550000 galaxies. This sample 
was divided at random into two subsamples with the same num- 
ber of objects generating a training and a validation set in order 
to calibrate the ANNz code used to the estimation of photometric 
redshifts. 

As a testing set for the photometric redshifts we used a random 
sample (hereafter Sz2) of 60,000 objects from MGS. We excluded 
from this sample galaxies belonging to the training set to avoid 
undesirable biases. 

F or the calibration of k-corrections we used the full MGS. Fol- 
lowing [MontMO^POTta&^rada] d2009l) . we have taken the magni- 
tude range where the number of galaxies per solid angle rises at a 
constant rate as a function of redshift in each SPSS band. There- 
fore, the apparent magnitude limits are set to ensure that the ef- 
fect of incompleteness is small in our sample. The adopted appar- 
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Table 1. Description of the samples used in this work. 



sample name 


number of objects 


Description 


Szl 


~ 550000 


Selected from SDSS-DR7 spectroscopic data. 

Used as training and validation sets in the computation of photometric redshifts 


Sz2 


~ 70000 


Selected from SDSS-DR7 MGS (excluding training set galaxies). 
Used as testing set for photometric redshifts. 


Sz3 


~ 82000 


Selected from SDSS-DR6 photometric data. 

Used to compare photometric redshift estimation with Oyaizu et al. (2008). 


Ski 


- 122000 


Selected from SDSS-DR7 MGS taking into account apparent magnitude and redshift limits (see text). 
Used for k-correction calibration. 


Sk2 


- 575000 


Selected from SDSS-DR7 photometric data with our photometric redshifts estimation. 
Used to compute k-correction at higher redshift and compare with literature. 




spec 



Figure 1. Left panel: z Bpcc vs z p hot relation for Sz2 sample (grey dots). The dashed line corresponds to the median value of the photometric redshift estimated 
at a given spectroscopic redshift interval. Error bars correspond to the 10 and 90 percentiles associated to the median. The solid line shows the one-to-one 
relation. Right panels: systematic (top right) and stochastic (bottom right) errors in the photometric redshift estimates. 



ent magnitude limits in each band are : u < 19.0, g < 17.91, 
r < 17.77, i < 17.24, and z < 16.97. For these limits we con- 
sider galaxies with z spec < 0.15 to assure completeness in the u 
and g bands, and z spec < 0.18 in the r, i and z bands. With these 
magnitude and redshift constraints, the sample used to calibrate the 
k-corrections (hereafter Ski) contains ~ 122000 galaxies. 

The photometric data set were extracted from the Galaxy ta- 
ble of the Cas Job 43 database. We restrict our analysis to pho- 
tometric objects with r < 21.5, since this magnitude limit as- 
sures good photometric quality and a reliable star-galaxy separation 
(Stoughton et al., 2002, Scranton et al., 2002). 

In order to contrast our estimates with those obtained by 
similar met hods, we compare o ur results with photometric red- 
shifts from lOvaizu et alj £2008). Taking into account these au- 



thors recommendation we use Photo zCC 2 estimates (obtaine d 
from colours and concentration indices, see lOvaizu et alj {2008)). 
It should be noted that Photoz2 table was not updated for SDSS- 
DR7, so we have constructed a sample from SDSS-DR6 photo- 
metric data selectingat random 82000 objects with Photo zCC2 
redshift information q For this sample we calculated photometric 
redshifts through the methods described in this work. Hereafter we 
will call this sample as Sz3. 

In order to compute k-correction at higher redshifts we have 
selected a random sample of ~ 575000 photometric SDSS-DR7 
data for which we have determined photometric redshifts and k- 
corrections (sample Sk2). For this sample we have also compared 



2 http://cas.sdss.org/dr7/ 



3 SDSS-DR6 photometric data was downloaded from Casjobs just as DR7 
photometric data including redshift information from Photoz2 table 
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the k-corrections computed in this work with different results ob- 
tained from the literature. 

Table Q] summarizes the main characteristics of the different 
samples used in this work 



3 PHOTOMETRIC REDSHIFTS 

Photometric redshift techniques use photometric parameters to per- 
form an estimation of galaxy redshift. This technique can be used 
to infer efficiently large numbers of galaxy distances, even for faint 
galaxies, for which spectroscopic measurements are prohibitive be- 
cause they would require large amounts of telescope time. 

There are different techniques to estimate photometric 
redshifts which can be classified into two groups. The first set of 
techniques makes use of a small number of model galaxy spectra 
derived from empirical or model-based spectral energy distribu- 
tions (SEDs). These methods estimate a galaxy redshift by finding 
an optimal combination of temp late spectra th a t reconstructs the 
observed galaxy colours (e.g., iBenitezI [2000; Bolzonell a et al. I 
120001 : ICsabaiet al.1l2003l) . The fact that these methods rely on a 
small number of template SEDs is their main disadvantage, in 
particular for galaxies at high redshifts, since a representative set of 
spectral templates applicable at all redshifts is not easy to obtain. 
The second group, called empirical methods (e.g.. lConnollv et all 
ll995tlBrunneretal.il 999). comprises techniques that need a large 
amount of prior redshift information, in general in the form of 
training set. This class of methods aims to derive a parametrization 
for the redshift as a function of photometric parameters. The form 
of this parametrization is obtained through the use of a suitably 
large and representative training set of galaxies for which we have 
both photometry and precisely known redshifts. In this case we 
can use combinations of galaxy photometric parameters, such 
as magnitudes in different photometric bands, galaxy colours, 
and concentration indices. The main drawback of empirical 
methods is that the training set should be representative of the sam- 
ple of galaxies for which we want to estimate photometric redshifts. 



3.1 Method 

We used Szl sample to compute photometr ic redshifts with the 
ANNz software package dCollister & Lahav 1 d2004l) ). which uses 
an Artificial Neural Network (ANN) to parametrize the relation be- 
tween redshift and photometric parameters. ANNz is based on a 
"multilayer perceptron"(MLP) algorithm, where the nodes are dis- 
posed in layers, and the nodes in a given layer are connected to all 
the nodes in adjacent layers. The ANN topology adopted in ANNz 
can be described as Ni n :N\:N2'----'-N out , where N% n and N ou t are, 
respectively, the number of input and output parameters, whereas 
Nj is the num ber of nodes in the i— th intermediate layer (e.g., 
Bishop1 [l99 lh . The first layer contains the inputs, which in our ap- 
plication are photometric parameters. The final layer contains the 
outputs, in this case the photometric redshift (z p hot). 

The free parameters of ANNz are the "weights" between the 
nodes, and are obtained by "training" the ANN with a training set 
consisting of galaxies with spectroscopic redshifts. The selected 
training set must be large enough and representative of the target 
population, to assure a reliable mapping of the input into the 
output. Also, it must contain the same set of input parameters than 
the target sample, for which we want to estimate the photometric 



redshifts. 

ANNz can be trained with different sets of input parameters 
in order to improve the photometric redshift accuracy. We have 
analysed several sets: i) magnitudes in the five SDSS bands, 
ii) colours and concentration index, iii) colours and Petrosian 
radii, and iv) SDSS-DR7 magnitudes in the five bands, plus 
concentration indices and Petrosian radii in g and r-bands. We 
found that the latter set is the best choice for redshift estimation. 
The use of concentration indices helps to break the degeneracies in 
the redshift-colour relation. This occurs due to the good correlation 
between the colour of a galaxy with its concentration indices and 
Petrosian radii. The resulting ANNz architecture adopted here is 
9 : 14 : 14 : 14 : 1. 



3.2 Results 

Our results of the ANNz training are shown in Figure Q] The left 
panel of Figure Q] shows the spectroscopic (z spcc ) vs. photometric 
(Zphot) redshift relation for the testing sample Sz2. As explained 
in Section 2, galaxies in the training set were excluded from the 
sample in order to avoid bias in this comparison. This figure also 
shows the one-to-one relation as a solid line, the median value of 
the photometric redshift estimated at a given spectroscopic red- 
shift interval (dashed line) and 10 and 90 percentiles (as error 
bars) of the scatter plot. The small dispersion seen in this figure 
corresponds to rms ~ .0227. For the MGS of the SDSS-EDR 
ICollister & Lahav I d20o4) obtained rms ~ 0.0229 and for a sam- 
ple comprising galaxies from SDSS-MGS, SDSS-L RG, CNOC2 
CFRS DEEP2 DEEP2, TKRS and 2SLAQ surveys. lOvaizu et all 
(2008) obtained rms ~ 0.03. The right panel of FigureQ]shows the 
systematic differences between jz p hot and z spcc , zu a s (higher sub- 
panel), and the rms dispersion a (lower sub-panel) as a function 
of redshift. Both, the bias and dispersion, shown little dependence 
with the spectroscopic redshift. 

F or c omparison with other e stimates, for LRGs. lAbdalla et al~l 
d2008h and lCollister et al. I d2007l) used the ANNZ code to compute 
a refined star/galaxy probability based on a range of photometric 
parameters. The photometric redshift rms deviation is 0.049 when 
averaged for all galaxies, and 0.040 for a bright sub-sample with 
i < 19.0 in the redshift range z p hot = 0.4 to 0.7. 

In Figure|2](top panel) we show the distribution of photomet- 
ric redshifts N(z) for Sz3 sample. In short dashed lines we display 
N(z) for our calculation, and in dotted lines Oyaizu et al. (2008) 
estimates. The solid curve corresponds to the expected distribution 
given by the theoretical calculation of Blanton et al. (2003), as- 
suming an universal luminosity function with parameters extracted 
from SDSS data. For z p hot < 0.1 the observed distributions are 
similar to the expected prediction. For photometric redshifts be- 
tween z p hot ~ 0.1 and z p hot ~ 0.35 our estimates are in better 
agreement with the theoretical curve, while Oyaizu et al. (2008) 
results show a ~ 20% deficit. We note that this range is very im- 
portant because SDSS-DR7 data has reliable completeness up to 
Zphot < 0.35, since beyond this redshift the galaxy density drops 
significantly. For z p i 10 t ~ 0.4 the observed distributions tend to 
overestimate the theoretical curve, an effect that could be associated 
to the Balmer break shifting between the q an d r filters, difficulting 
the redshift estimates dBudavari et al. I d200lh). The bo t tom p anel of 
Figure [2] show the comparison between I Oyaizu et "all d2008h •^phot 
and our estimates. 
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Figure 2. Top panel: distribution of photometric redshifts N(z) for our calculation (short dashed) and Oyaizu et al. 2008 (dotted) estimates. The solid curve 
corresponds to the expected distribution given by the theoretical calculation of Blanton et al. 2003. Bottom panel: the 2 p hot redshift relation for Oyaizu et al. 
(2008) estimates and our estimates. The dashed line corresponds to the median value and error bars correspond to the 10 and 90 percentiles associated to the 
median. 




Figure 3. k-correction estimated by different authors. The solid line corresponds to the median value, dashed line show the one to one relation and error bars 
correspond to the 10 and 90 percentiles associated to the median. The large scatter between different methods can be appreciated. 
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Figure 4. Median K-correction estimates for the 5 SDSS photometric bands u, g, r, i and z. The lines represent the fit to the median for 5 cuts in reference- 
colour index (g — r). The error bars corresponds 25 and 75 percentiles from the median. The dotted lines in k g and k r represent the initial k values for the 
iteration (see text). 



In Figure [2] (bottom panel) shows the z p hot redshift relation 
for Oy aizu et al. (2008) esti mates and our estimates. 

In lO'MilletalJ ( l2008l) it is estimated the maximum redshift 
(zmax) beyond which compact galaxies are confused with the im- 
age PSF. For SDSS, a galaxy with M r = —23 has a z max of 0.63 
(see lO'Mill eTai] J2008T) Table 1). For this reason, we recommend 
using galaxies with z p hot ^ 0.6. 



4 K CORRECTIONS 

k-corrections allow to transform the observed magnitudes at a red- 
shift z, into standard luminosities at some reference frame. This 
correction depends on the filter that was used for the observations, 
the rest-fra me standard, t he sha pe of the SED of the galaxy, and the 
redshift (see lHogg et "all ( 1200 lh ). 

In the previous sections we have used Artificial Neural Net- 
work as a tool to perform non-linear fitting to obtain photometric 



redshifts. An analysis based on a neural network approximation is 
the natural trend to obtain k-corrections. However these techniques 
needs a large and representative amount of prior information in the 
form of a training set. 

The problem with the former is the lack of confidence on indi- 
vidual k-corrections estimates. This is shown in Figure [3] where 
we compare r-band k-corrections estimated by different authors for 
6113 galaxies, selected randoml y from MGS, in t he narrow red- 
shift range 0.14<z<0.1 6. ForlBlanton&Ro weis (2007) we use 
k-correct_v4 . 2, for Chilinaarian et al. I 112010]) we employed 
the on-line available codeK for lWestra et al.ld201fj|) we used the on- 
line empirical k-correction calculator and using models and (g- 
r) parameter and we also compa re with k-corrections downloaded 
from the SDSS Photoz table ( Aba zaiian et al. I d2009l) ). 



4 http://kcor.sai.msu.ru/getthecode/ 

5 http://tdc-www.cfa.harvard.edu/instruments/hectospec/progs/EOK/ 
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Figure 5. A and B coefficients for the K correction fits of the form k = Az + B where A = a a {g — r) + b^ and B = ag(g — r) + bg 



It is clear that individual estimates have a large uncertainty as 
indicated by the large scatter between different methods, whereas 
the mean trends can be trusted. 

Since the deve lopment of the public code k-correc{f| 
( Blant onet al. 1 12003)) many authors have used this method to cal- 
culate k-corrections. This code is based on a mathematical algo- 
rithm, namely a non-negative matrix factorization, which creates 
model based template sets. The set of templates is reduced to a ba- 
sis of the five SDSS pass bands, which are used to interpret the SED 
of the galaxy in terms of stellar populations. Then, linear combina- 
tions of these templates are used to fit spectral energy distributions 
to broadband photometric observations for each galaxy and so k- 
corrections are obtained. 

The disadvantage of this method is the use of a limited number 
of spectra that decreases with distance. The intrinsic colour of a 
galaxy is related to its SED, therefore this property can be used 
to compute k-corrections at high redshifts, where the Blanton's k- 
correction technique is more uncertain. 

In this paper we use the public available code described in 
iBlanton & Roweisl j2007l) (k-correct_v4 . 2) as a tool for cali- 
brate our k-corrections. 

4.1 Model for k-corrections 

The use of galaxy colours to describe galaxy populations presents 
the advantage that colours are easily quantifiable rather than mor- 
phological types. In this work we analyse a relation between k- 
corrections and the reference-frame galaxy (g — r) colour index. 

As described in Section 2, Ski sample comprises MGS up 
to z ~ 0.2. Since the motivation of this work is to derive k- 
corrections at higher redshifts, we have calculated k-corrections 
for the Ski sample using the k-correct_v4 . 2 code, analysed 

6 http://howdy.physics.nyu.edu/index.php/Kcorrect 



their dependence on the galaxy reference frame (g — r) colour 
and extrapolate the results obtained to higher redshifts. Follow- 
ing IBlanton & Roweis ] d2007l) , k-corrections were calculated on 
the five SDSS photometric bands shifted to z apec = 0.1. 

We divide the spectroscopic data into five different sub- 
samples (d) according to the reference-frame (g — r) colour: 

Ci:(fl-r) <0.5, 
C 2 :0.5< (g-r) < 0.65, 
C 3 : 0.65 < (g - r) < 0.8, 
C 4 : 0.8 <{g-r) < 0.95, 
C 5 :(g-r) > 0.95. 

Then we calculate median values of k-corrections per redshift 
interval for the different C\ sub-samples, for the five SDSS bands. 
The results are shown in Figure [4] where it can be appreciated the 
smooth dependence of the median k-correction values as a function 
of redshift once the reference galaxy (g — r) colour is taken into 
account. We consider galaxies with z ^ 0.1 and we adopt a linear 
relation k — Az + B for each d sub-sample. The resulting fits 
are also displayed in Figure [5] where the error bars correspond to 
25 and 75 percentiles of the median values. 

By inspection of this figure it can be seen that our linear model 
is a good representation of the relation between k-correction and 
redshift. 

The values of the parameters A and B are obtained by fitting 
a straight line through \ 2 minimization for each d sub-sample. 
In Figure [5] we plot the derived values of A (upper panel) and B 
(bottom panel) as a function of the mean (g — r) values for each 
d sub-sample; the error bars correspond to the square root of the 
variances in the estimates of A and B from x 2 minimization. From 
this figure it can been seen that a linear relation also gives a fair 
representation of the trends observed. Therefore, we adopt a model 
where A = cia(<? — r) + 6.4 and B — ag(g — r) + bg- The final 
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Figure 6. k-correction obtained for a random sample of SDSS-DR7 galaxies. Our estimates are sh own in grey and in dark-grey we provide the results using 
Blantons 2003 code. Blue solid li nes are k-corrections obtai ned for early-type galaxies presented in lRocheetalJi2009». red solid lines k-correction computed 
using "K-correction calculator" JChilingarian et al. 1 120101) ). magenta solid lines corresponds to Abazai ian et al. 1 120091) estimates and in green solid lines 
IWestra et alJfcOlOh k-corrections. 



calibration of Blanton's k-corrections is given as 

kj = [aAj(g -r) + b A j]z + [a B] {g - r) + b B j] (1) 

where j represents the different ugriz bands. 

In Table [2] we present the a^j, b A j, a B j and b B j values and 
the errors obtained by % minimization. 



4.2 Estimation of k-corrections 

In order to compute k-corrections using the model described in the 
previous Section, it is necessary to obtain reference-frames (g — r) 
colours for each galaxy, which in turn depends on k-corrections. 

To avoid this problem we have adopted an iterative procedure 
as follows: the iteration starts with an initial value for k g and k r 
according to the value of the concentration index of the galaxy in 
the r band (c r ). This quantity does not depend on k-correction and is 



a suitable indicator of galaxy morphological t ype bimodality: early - 
type galaxies hav e higher c r than later types i Strateva et al. ifeOOlh. 



Kauff mann et alj d2003d) . iKauffmann et al.l J2003 j) , iMateus et all 
(2003)). 

According to Strateva et alj ( l200ll) . c r > 2.55 values cor- 
respond to early-type galaxies, and late-type galaxies have c r < 
2.55. This bimodality can also be seen in the colour distribution 
of galaxies. In particular, the distribution of (g — r) colours has 
two well-defined peaks: one at (g — r) — 0.60 and other at 
(g — r) = 0.95, corresponding to the late and early-type compo- 
nents, respectively. Taken these properties into account, we assign 
(g - r) initial = 0.60 if c r < 2.55 and (g - r) ini uai = 0.95 if 
c r > 2.55 so that the initial k g and k r values are: 

kl,initial = [a-Al(g — r)i n itial +bAl]z+ [asi(g — r)i n itial +&Bi], 

where I refers to the g and r bands (see the dotted lines Figure|4j 
Once these initial values are fixed, we iterate in equationQ]to 



Photometric redshifts and K-corrections for SDSS-DR7 9 




Table 2. Parameters of the kj = [a,Aj (g — r) + b \j\z + [asj (g — r) + bgj] relation obtained by \ 2 minimization 



Band 




<Ta A 


b A 




a B 


<?a B 


b B 




u 


2.956 


0.070 


-0.100 


0.034 


-0.299 


0.009 


-0.095 


0.004 


9 


3.070 


0.165 


0.727 


0.117 


-0.313 


0.021 


-0.173 


0.015 


r 


1.771 


0.032 


-0.529 


0.023 


-0.179 


0.005 


-0.048 


0.003 


i 


0.538 


0.085 


-0.075 


0.079 


-0.027 


0.013 


-0.120 


0.012 


z 


0.610 


0.045 


-0.064 


0.034 


-0.061 


0.007 


-0.106 


0.005 



obtain k g and k r for each galaxy. After each iteration, we check 
that the colours are within an acceptable range, < (g — r) < 1.8 
in order to avoid either too red or too blue colours. In these cases 
the iteration starts again with a new (g — r) ini t ia i value selected at 
random from a Gaussian distribution (within lcr) that fits either the 
blue or the red peak of the reference frame colour distribution. 

The iterations stop when the difference in both k g and k r be- 
tween two consecutive steps is less than 0.001. This procedure con- 
verges in less than 15 iterations. 

From the finally obtained k g and k r values we calculate the 
reference frame (g — r) colours, which allows to compute k- 
corrections in the other bands using equation Q] 

We notice that there is a small percentage (less than 0.4%) 
of galaxies for which our algorithm does not converge. Howerver, 
the main sources for this lack of convergence are large magnitude 
uncertainties and unreliable observed colours((g — r) b s > 3). 
The last colours constraints are helpful on remo ving stars wit h 
unusual colours, wi t hout di scarding real ga l axies jLopesI J2007I) . 
Padmanabh an et alj J2005h ). ICollister et all d2007h point out that 
the stellar contamination may still be present. 



4.3 Results 

We have compared the k-corrections obtained in this work with 
the results from k-correct_v4 . 2 using the Sk2 sample. As ex- 
plained in Section 2, we have restricted our analysis to galaxies 
withr < 21.5. 

In Figure [6] we plot both estimations against redshift for the 
five SDSS-DR7 bands. In dark-grey we show k-correct_v4 . 2 
results whereas light-grey points are our estimates. It can bee seen 
in all cases, that ou r k-correction shows the s ame trend but with 
a lower spread than lBlanton & Roweis 1 (l2007h results. This lower 
spread is particularly noticeable in the i and z bands. However , 
our estimations are in agreement with IChilingarian et al. I (2010) 
k-corrections in these bands. In Figure [6] we also co mpare our k- 
correc tions with those obtained for E/S0 galaxies by Roc he etHll 
(2009). We find that the mean k-correction for the early type 
galaxies lies within our estimates for red galaxies (Cs sample), 
where k-correction values are higher. We also perform a compari- 
son with k-correctio n s obta ined using the "K-correction calculator" 
( Chilingarian et al. the "on-line empirical k-correction 

calcul ator" dWestra et alj|201fj|) and k-corrections from Photo z 
table dAbazaiian et al. 1 120091) . We compute the k-corrections in 
all bands at extreme red and blue colours, finding that these esti- 
mates are within our range of calibrated values. We notice a differ- 
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Figure 8. Luminosity distribution for galaxies brighter than M r — 5log(hi^) = —21.5 and for 4 different redshifts cuts. 
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Figur e 9. Left panel: Lumin osity function in g and r bands. The dashed lines show the Schecter fit for MGS of iBlanton et al. I J2003l) and 
Montero-Do rta & Pradal J2009I) ( green and blue respectively). The green solid line is the result from the l/Vrnax method for our photometric s ample . 
The red solid lines are the result from the STY method from two different redshift bins (0.1 < z < 0.35 and 0.35 < z < 0.55) from lZucca et alj J20O9I) . 
The error bars (blue colours) represent lcr uncertainty calculated using a bootstrapping technique. Right panels: V/Vmax test. The error bars represent lcr 
deviation from the median. The horizontal dotted lines at V/Vmax = 0.5 correspond to the g and r band completeness. 



ence in the k u for blue g alaxies at higher redshift with respect to 
Chilingaria n et al.l ( 120101) . We argue that this can be originated in 
the fact that these authors calibrate k u with the u — r colour and 
that the it filter has a natural red leak that cases abnormal u — r 
colours. This effect can propagate to the k-corrections, generating 
higher values, particularly at large redshifts. 



In Figure [7] we plot the distributions of k g and k r in four 
different redshift intervals. Our estimates are shown as solid 
lines, and the dashed h istogram is k-correct_v4 . 2 code 
(Blanton & Roweis 120071) results. For galaxies in the redshift range 
0.1 < Zphot < 0.3, both distributions are similar. As redshift in- 
creases, our distributions remain always unimodal, with a well de- 
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Figure 10. Derived rest frame (g — r) colour distribution for galaxies brighter than M, — 5iog(/i7s) = —21.5 and for 4 different redshifts cuts in Sk2 sample. 
The dashed histogram corresponds to Blanton's 2003 results and in solid lines our estimates. Relative excess of red galaxies in Blantons results is clearly seen 
in comparison to our colour estimates 



fined mean. On the other hand, Blanton's k-correction distribution 
shows 3 different maxima, one centered approximately in our dis- 
tribution, and the others in the extremes. This behavior is probably 
indicative of template mismatch. At redshifts 0.1 < 2 p hot < 0.6, 
both k r distributions approximately match each other, whereas g- 
band k-corrections from k-correct_v4 .2 distribution is lop- 
sided to negative values. 

4.4 Luminosities and colours 

We have estimated the absolute magnitude of galaxies in the Sk2 
sample in the g and r bands using our k-corrections. The absolute 
magnitude of a galaxy in a given band is related to its apparent 
magnitude by: 

Mi = l- 5log 10 (DL(z)) - 25 - k t 

where I refers to the g and r bands, and DL(z) is the luminosity 
distance (which depends on the cosmological parameters adopted). 
In Figure [8] we show the distribution of k-corrected absolute mag- 
nitudes in four different increasing redshift intervals (key in Fig- 
ure). The shape of these distributions is very similar in both g and 
r SDSS bands: a bell-shaped distribution skewed to fainter mag- 
nitudes. Notice that the mean of the distributions move towards 
fainter magnitudes as the redshift range decreases, while the bright 
tail of the distribution remains approximately fixed at M g ~ —22.5 
and M r ~ -23. 

For Sk2 photometric sample we computed the luminosity 
function in the g and r bands using the 1 /Vma x method consi d- 
ering the incompleteness with a V/Vmax test ([Schmidt ( 1968)). 
This methods takes into account the volume of the survey enclosed 
by the galaxy redshift and the difference between the maximum and 
minimum volumes within which it can be observed. 

In Figure [9] we show the luminosity function and V/Vmax 



test for the redshift range 0.1 < z P hot < 0-6. The da shed 
lines show the Schechter fit for MGS of lBlanton et al. I d2003l) and 
iMontero-Dorta & Prada I d2009l) . The green solid line is the result- 
ing Schechter fit for Sk2 sample. We compare d our result in g ban d 
(4686 A) with two different redshift bins from lZucca et alj {2009). 
These authors have studied the evolution in the B band (4459. 7 A) 
luminosity function to redshift z ~ 1 in the zCOSMOS from the 
STY method. 

Our best Schechter fit corresponds to M* ~ —20.469 ± 
0.0535 , $* ~ 0.0224 ± 0.0089 and a ~ -1.065 ± 0.0286 for g 
band, and M* ~ -20.821 ± 0.0966, $* ~ 0.030 ± 0.0029 and 
a ~ — 0.78 ± 0.0298 for r band. Where the errors had been ob- 
tained through bootstrapping. We notice th at the luminosity func- 
tions at 2 ~ 0.5 from our work and that of IZucca et all 12009) arc 
consistent taking into a ccount the (g — B) values of typical Sbc 
galaxies at this redshift (Fukug ita et aljl99a) . 

In Figure [10] we compare the (g — r) colour distribution of 
galaxies brighter than M r = —21.5 in four redshift intervals. 
The dashed histogram corresponds to k-correct_v4 . 2 code 
results, and in solid lines we show our estimates. For galaxies 
with 0.1 < Zphot < 0.3 both distributions are similar, exhibit- 
ing a prominent red peak corresponding to early-type galaxies. As 
redshift increases, the distributions are shifted bluewards, and we 
notice in the redshift range 0.45 < Zp hot < 0.6 an excess of 
red galaxies in lBlanton & Roweis 1 d2007h results in comparison to 
ours. This trend can also be seen in the full redshift range (Figure 
[To] bottom right panel), where the galaxy bimodality can also be 
appreciated. 

In Figures QT] and Q2] we show the colour-colour diagrams 
(g — r) vs (r — i) and (g — r) vs (i — z) for galaxies brighter than 
M r = —21.5 in four different redshift intervals. The contours en- 
close 50%, 65% 85%, and 90% of the galaxies in Sk2 sample. In all 
cases the shaded region represents the corresponding colour-colour 
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Figure 11. (g ~ r) vs (r — i). colour-colour diagram for galaxies brighter than M r = —21.5 in four different redshift intervals for Sk2 sample. The contours 
enclose 50%, 65% 85%, and 90% of the photometric sample. The shaded region represents the colour-colour diagram for 90% of MGS. 



diagram for 90% of MGS, which have spectroscopic redshifts. The 
left upper panel (0.1 < Zphot < 0.2) shows that contours cal- 
culated for the photometric sample matches the shaded region. As 
the redshift increases, the contours in the colour-colour diagrams 
expand with respect to the spectroscopic sample, but both have ap- 
proximately the same center. The expansion of the contours for the 
photometric results could be due to uncertainties in the determi- 
nation of photometric redshifts and k-corrections, as well as they 
may reflect galaxy evolution. The colours show a decreasing trend 
for the red populati on as well as a constant shift to bluer colours 
dO'Mill et all 00081) ). This is particularly important for the redshift 
range 0.45 < z phot < 0.6. 



5 CONCLUSIONS 

In this work we present a new set of photometric redshift (2 p hot) 
and k-correction estimations for the SDSS-DR7 photometric cat- 
alogue available on the World Wide Web. In order to calculate 



Zphot, artificial neural networks were applied using the public code 
ANNz. The improvements in the SDSS-DR7 photometric redshift 
estimation are: 

1) We added the concentration index and the Petrossian radii in g 
and r bands to the usual five magnitudes used in previous similar 
works. These additional inputs improve z p hot estimations because 
the concentration index provides information regarding the slope of 
the galaxy brightness profile, helping us to break the degeneracies 
in the redshift-colour relation due to morphology. The Petrossian 
radius is a robust measure of how shallow the brightness profile is 
and contain information about the angular size, that is related to the 
distance. 

2) The choice of different galaxy samples for the training set (MGS, 
LRG and AGN sample) provides a wide sampling of different 
galaxy types at various redshifts, allowing to improve z p hot esti- 
mates. 



Our Zphot estimates have a rms ~ 0.0227, and the resulting 
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Figure 12. (g — r) vs (i — z) colour-colour diagram for galaxies brighter than M r = —21.5 in four different redshift intervals for Sk2 sample. The contours 
enclose 50%, 65% 85%, and 90% of the photometric sample. The shaded region represents the colour-colour diagram for 90% of MGS. 



galaxy distribution shows a good agreement with the theoretical 
distribution derived from the SDSS galaxy luminosity functions. 

We have used the k-correct_v4 . 2 code 
(Blanton & Roweis for the MGS and we have per- 

formed a linear fit between reference frame (g — r) colour and 
redshift, extrapolating this relation at high redshifts. We propose 
an iterative procedure to estimate k-corrections from the observed 
photometry in the g and r bands. Using initial values that depend 
on the concentration index and the observed colour, we obtain the 
k-correction for the other bands. Our results show that the use of 
this simple linear relation between the reference frame (g — r) 
colour and redshift is as accurate as those obtained in previous 
work. A clear plus of our approach is the low computational time. 

Our k-correction estimations do not use templates, avoiding 
statistical errors in the lack of homogeneity in spectral features, and 
minimizing systematical errors caused by an assumed spectral en- 
ergy distribution (SED). This can be noticed in the smooth behavior 
of the distribution of k-corrections, even for intermediate redshifts. 



The analysis of the distribution of k-corrected absolute mag- 
nitudes show that the shape of these distributions has a bell-shaped 
skewed to fainter magnitudes and the mean of the distributions 
move towards fainter magnitudes as the redshift range decreases. 

We have computed the luminosity function in g and r bands 
trhough 1 /Vmax me thod taking into account the incompleteness 
with a V I Vmax test dSchmidtl jl968hl. We notice that the curves 
derived from this work and that of lZucca et al.|j2009l) are consistent 
taking i nto account the (g — B) value of a typical Sbc galaxy at 
2 ~ 0.5 lFukugita et alj dl995h . 

The analysis of the (g — r) colour distribution for galaxies 
brighter than M r = —21.5 shows that the galaxies in our samples 
are shifted bluewards as redshift increases. This trend leads to the 
emergence of bimodality in the full redshift range. From the colour- 
colour diagrams, we can conclude that the behavior of colours at 
low redshift is in good agreement with the trends of the spectro- 
scopic sample. As redshift increases we see a broadening of the 
contours and an increase in the blue galaxy population. However, 
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the distribution of galaxies in the colour-colour diagram remains 
centered with respect to the spectroscopic data. 
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